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Abstract 

In the literature, there are five distinct, fragmented sets of analytic predictions 
for the scaling behaviour at the phase transition in the random-site Ising model in 
four dimensions. Here, the scaling relations for logarithmic corrections are used to 
complete the scaling pictures for each set. A numerical approach is then used to con- 
firm the leading scaling picture coming from these predictions and to discriminate 
between them at the level of logarithmic corrections. 



1 Introduction 



One of the major achievements of statistical physics is the fundamental explanation of 
critical behaviour at continuous phase transitions through Wilson's renormalization group 
(RG) approach. While this has mostly provided a satisfying picture for over thirty years, 
certain types of phase transitions have resisted full treatment. Such stubborn cases, which 
have been the subject of conflicting proposals and analyses, include systems with in-built 
disorder. 

The Ising model with uncorrelated, quenched random-site or random-bond disorder 
is a classic example of such systems and has been controversial in both two and four 
dimensions. In these dimensions, the leading exponent a which characterises the specific- 
heat critical behaviour vanishes and no Harris prediction for the consequences of quenched 
disorder can be made [1] . The Harris criterion indicates if the specific heat of the pure 
system diverges - i.e., if a > there - then the critical exponents may change as random 
quenched disorder is added to a system. If a < in the pure system, then this type 
of disorder does not alter critical behaviour and the critical exponents are unchanged in 
the random system. In the two-dimensional case, the controversy concerns the strong 
universality hypothesis which maintains that the leading critical exponents remain the 
same as in the pure case and the weak universality hypothesis, which favours dilution- 
dependent leading critical exponents (see [2] and references therein). 

Since ci = 4 marks the upper critical dimensionality of the model, the leading critical 
exponents there must be given by mean field theory and there is no weak universality 
hypothesis. However, unusual corrections to scaling characterise this model, and the 
precise nature of these corrections has been debated. This debate motivates the work 
presented herein: methods similar to those employed in [2], namely a high-statistics 
Monte Carlo approach coupled with finite-size scaling (FSS), are used to progress our 
understanding of the four- dimensional version of the random-site Ising model (RSIM). 

While not directly experimentally accessable, the four- dimensional RSIM is of interest 
for the following reasons: (i) it is closely related to the experimentally important dipolar 
Ising systems in three dimensions, (ii) it is an important testing ground for the widespread 
applicability of the RG, (iii) it presents unusual corrections to scaling, (iv) in high energy 
physics, the establishment of a non-trivial Higgs sector [3] for the standard model requires 
a non-Gaussian fixed point and a new universality class which may, in principle, result 
from site dilution and (v) it is the subject of at least five analytical papers which differ 
in the detail of the scaling behaviour at the phase transition. 

In Section [21 the current status of the RSIM in four dimensions is reviewed, paying 
particular attention to previous detailed analytical predictions in the literature for its 
scaling behaviour. The scaling relations recently presented in [H [5] are then used to 
construct full scaling descriptions based on earlier partial theories [HI [TJ [H [9l [10]. In 
Section [3] the theoretical finite-size scaling behaviour of the model is presented. The 
details of the extensive numerical simulations are given in Section [H the outcomes of 
which are analysed in Section [51 Our conclusions are summarised in Section [61 
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2 Scaling in the RSIM in Four Dimensions 



The consensus in the hterature is that the following structure characterises the scaling 
behaviour of the specific heat, the susceptibility and the correlation length at the second- 
order phase transition in the RSIM in four dimensions (up to higher-order correction to 
scaling terms) [HI [3 [H El [TOl [11] : 



Coo(t) 



Xooit) ~ 



eoc(t) ~ 

Here, the subscript indicates the size of the system, the reduced temperature t = {T — 
Tc)/Tc marks the distance of the temperature T from its critical value Tc and A and B > 
are constants. The correlation function at criticality decays as [TJ [9] 

Gooix) = x-^'^-^+'^^\lnx\\ (2.4) 

where x measures distance across the lattice, the dimensionality of which is d. The cor- 
relation length for a system of finite linear extent L also exhibits a logarithmic correction 
and is of the form 

eL(t = 0) ~L(lnL)^ (2.5) 

The leading power-law behaviour is believed to be mean field because the fixed point is 
expected to be Gaussian and 

a = 0, (3=^, 7 = 1, 5 = 3, u = ^, r/ = 0, A = ^. (2.6) 

Here, f3 and 5 are, in standard notation, the critical exponents for the magnetization out 
of field and in field, respectively while A is the gap exponent characterising the Yang-Lee 
edge. There is no dispute in the literature regarding these leading exponents, some of 
which will be re-verified in this work. Nor is there any dispute regarding the details of 
the unusual exponential correction terms in f l2.1l) -( |231) . However there are at least five 
different sets of predictions for the exponents of the logarithmic terms, which differ from 
their counterparts in the pure model, and a principle aim of this work is to investigate 
these predictions numerically. 

Aharony used a two-loop renormalization-group analysis to derive the unusual expo- 
nential terms in fl2.1l) - fl2.3p . and also found [B] 

d = ^, 7 = 0, u = 0. (2.7) 



A-5|t|-"exp I -2W— |ln|t|| | | In |t| 



\tr'exp\\l^\\n\t\\\\\ii\t\\\ 



ir'^expiy^llnltlllllnltir 



(2.1) 
(2.2) 
(2.3) 
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In [7], Shalaev pointed out that Aharony's results needed to be refined and, by determining 
the beta function to three loops, gave predictions for the specific heat and the susceptibility 
which differ from those in [6] in the slowly varying multiplicative logarithmic factors: 

d = 1.2368, 7 = -0.3684, r) = 0094 . (2.8) 

Jug studied the a = line of ra-component spin models in (n, d) space where d is the 
system's dimensionality, and thereby worked out the logarithmic corrections for the d = 4 
n- vector model [8]. For the case at hand {n = 1), he obtained 

a = 1/2, 7 = 1/212 ^ 0.0047 . (2.9) 

In [9], Geldart and De'Bell confirmed that to obtain the correct powers of | In |t|| the beta 
function has to be calculated to three loops, but the results of [9] differ from those of [7] 
in the powers of the logarithms which appear in the specific heat and in the correlation 
function: 

a ^ 1.2463, 7 ^ -0.3684, V= ^ = 0.0047. (2.10) 

Finally Ballesteros et al. [10] extended and corrected Aharony's computation to give the 
correction exponents: 

a = -, 7 = — ^ 0.0094, = 0, q = -. (2.11) 
2' ' 106 ' ^8 ^ ' 

So the detailed analytic scaling predictions of at least five groups of authors clash, 
and a number of questions arise: (i) are the predictions from within each author set 
self-consistent? (ii) what are the full set of predictions (i.e., extended to all observables) 
coming from each set? (iii) can a simulational approach yield numerical support for the 
shift in the correction terms from their counterparts in the pure model, and, further, 
(iv) can such a computational approach support one or other of these five different sets 
of analytic predictions? Here the scaling relation for logarithmic corrections developed 
in [U [S] are used to accomplish (ii) and it is shown that the answers to questions (i) and 
(iii) and to some extent (iv) are in the affirmative. In particular, numerical support for 
the broad scenarios presented in [6l [H [10] is presented. 

Modification of the self-consistent scaling theory for logarithmic corrections of [H [5] , 
to incorporate the exponential terms, leads to the following forms for the behaviour of 
the magnetization in the 4D RSIM: 



moo(t) = t'^exp |^-i^^|ln|t||j|lnt|^ (2.12) 

m^{h) = h^\nh\^ . (2.13) 
From Eq.(15) of [4j, we also write for the scaling of the Yang-Lee edge 



rYL(t) -t'^exp I -^Jl|ln|t|| l|lnt|^. (2.14) 
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Table 1: Theoretical predictions for the exponents of the logarithmic corrections to scaling 
for the pure Ising model in four dimensions and for its random-site counterpart. The 
latter exponents are listed in boldface if they come directly from the cited literature. 
The remaining values are extended from those of the literature using the scaling relations 



Log 


Pure model 


Aharony [6j 


Shalaev 12] 


Jug [8j 


Geldart 


Ballesteros 


exp 










& De'Bell [9] 


et al [10] 


a 


1/3 


0.5 


1.237 


0.5 


1.246 


0.5 




1/3 


0.25 


0.434 


0.252 


0.439 


0.255 


1 


1/3 





-0.368 


0.005 


-0.368 


0.009 


5 


1/3 


0.167 


0.167 


0.170 


0.170 


0.173 


V 


1/6 





-0.189 




-0.187 





V 








0.009 




0.005 


0.009 


q 


1/4 


0.125 


0.120 




0.125 


0.125 


A 





0.25 


0.803 


0.248 


0.807 


0.245 



The scaling relations for logarithmic corrections in this 4D model are [H [5] 



a = dq — do , (2.15) 

2/3-7 = dq-du, (2.16) 

/3(5-l) = 66-^, (2.17) 

V = l-0{2-r]), (2.18) 

A = /5-7. (2.19) 



These scaling relations are now used to generate a complete scaling picture from the 
fragments available in the literature O [3, [H [9l [10]. This complete picture is given in 
Table [1], where the exponents of the logarithmic correction terms are listed. Values for 
the exponents in boldface come directly from the reference concerned and the remaining 
values are consequences of the scaling relations ( 12.15p -( l2.19p . Each of the five papers P, 
[3, El El [To] is self-consistent in that the exponents given within do not violate logarithmic 
scaling relations. However, there are clear discrepancies between each of the five papers. 

The presence of the special exponential corrections has recently been verified by Hell- 
mund and Janke in the case of the susceptibility [11]. These exponential terms mask the 
purely logarithmic corrections, so in order to detect and measure the latter one needs to 
cancel the former. Certain combinations of thermodynamic functions achieve this, but it 
turns out that FSS does this also. FSS therefore offers an ideal method to determine the 
exponents of the logarithmic corrections numerically [2]. 

^ The relation (|2.15p is modified to read a — 1 + dq — diy when a = and when the impact angle of 
Fisher zeros onto the real axis is any value other than 7r/4, which is not expected to be the case in this 
4D model [5]. 
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3 Finite-Size Scaling 

Fixing the ratio of ^oo(i^) in ( 12. 3p and ^l(O) in ( 12. 5p to x, one has 



2 V 53 

Taking logarithms of both sides, one obtains 



t-'exp I ij|-|ln|t|| ||ln|t||'' = xL(lnL)^ (3.1) 



ln|t||^-lnL, (3.2) 

V 



which re-inserted into fl3.ll) gives 



t ^ L-tduD^expli^lilnL lH-oUL^ (3.3) 



- i-^(>nir*exp(^H,„ij{i + o(_i^)}. (3.4) 

having used the mean-field value (12.61) for the leading exponent v and the logarithmic 
scaling relation (12.151) . If a = 1/2, this recovers a result in [10] for the FSS of the 
pseudocritical point. 

Inserting (13.41) into (12.31) recovers (12.51) . as it should. The FSS of the remaining func- 
tions are determined by inserting (I33D into (O) to (EJD and fl212D to fl2AiD . One finds 



Cl{0) ^A- B'L^ exp ( - ( 2 + — ) J— In L \ (In L)"+^«-'') , (3.5) 



2vJ V 53i/ 

where B' (x B is & positive constant [HI [71 [H [U [10] . Inserting the mean-field values a = 0, 
V = 1/2, one obtains the simpler form 



Cl{0) ^A-B' exp 1^-2^ 1^ In L j (In Lf . (3.6) 

Similarly, the FSS for the susceptibility is 

Xl(0) ~ L^\\nL[<-'-^^-^^ = L^llnLlS (3.7) 

where 

C = 7-2(z>-g) = ^d + 7. (3.8) 

That for the Yang-Lee edge is 

n{L) ~ L-^\\nL\^+^^'''^^ = L-^\\nL\P, (3.9) 
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Table 2: The exponents of the multiphcative logarithmic corrections to FSS for the mag- 
netic susceptibility and for the Lee- Yang zeros coming from the literature and compared 
to their equivalents in the pure case. The FSS exponents are ( for the susceptibihty and 
p for the Yang-Lee edge. 



Exponent 


Pure 
model 


Aharony 

[6] 


Shalaev 


Jug 

i8j 


Geldart & 
De'Bell [9] 


Ballesteros 
et al pO] 


Susceptibility ( 
Lee- Yang zeros p 


1/2 
-1/4 


0.25 
-0.125 


0.25 
-0.125 


0.255 
-0.127 


0.255 
-0.127 


0.259 
-0.130 



where 

p = A + 3(z>-g) = -^d-^7- (3.10) 

Each of these also have sub-leading scaling corrections of strength (9(l/\/lnL) times the 
lead behaviour. One notes, however, that the unusual exponential terms, which swamp 
the logarithmic corrections in the thermal scaling formulae (12.21) and fl2.14p . drop out of 
their FSS counterparts (13.71) and (13.91) . These are therefore ideal quantities to study the 
logarithmic corrections. The theoretical analytical predictions of each of the five sources 
in the literature are now used to construct five possible FSS scenarios for the specific 
heat, the susceptibility and the Lee- Yang zeros. While Jug did not calculate the critical 
correlator or correlation length in 4D, the FSS picture corresponding to [8] can still be 
constructed through the scaling relations for logarithmic corrections. The FSS scenarios 
are listed in Table O 

The remainder of this paper is concerned with Tables [U and El The primary objective 
is to verify that the exponents for the logarithmic-correction terms in the RSIM are indeed 
different to those of the pure model. Once this is established, one would like to determine 
which of the five sets of analytical predictions are supported numerically. From Table 2, 
it is clear that present-day numerics cannot be sensitive enough to distinguish between all 
five scenarios for the susceptibility or individual zeros. However, there are clear differences 
between the predictions coming from [6l[8l[T0] and [7l[9] for the specific heat (Table 1) and 
it will turn out that the numerical data is indeed sensitive enough to favour the former 
over the latter. 

4 Simulation of the RSIM in Four Dimensions at Var- 
ious Dilution Levels 

We have performed extensive simulations of the model for linear lattice sizes from L = 8 
to L = 48 at dilutions p = 1, p = 0.8 and p = 0.5. In each case, we have employed a 
Wolff single-cluster algorithm [13] to update the spin variables using periodic boundary 
conditions. Thermalization tests including the comparison between cold (all spins up) 
and hot (all spins random) starts have been carried out. We found that the plateau 
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Figure 1: (Color online) Averaged behaviour of the susceptibility with the Monte Carlo 
time for 20 samples at L = 32 and p = 0.800. After every Monte Carlo sweep (Wolff 
update) measurements were performed. The plateau is reached more easily starting from 
a cold configuration (triangles). 



for the susceptibility is quickly reached by starting from cold configurations, see Fig. [H 
Indeed, the results for the susceptibilities from hot and cold starts are fully compatible 
(and are less than two standard deviations away from each other, even at the level of 
logarithmic corrections). The information about the numerical details is given in Table [31 
We have taken 1000 disorder realizations in all the cases except for L = 48, where only 
800 samples were used. We estimate that the total simulation time has been equivalent 
to 20 years of a single node of a Pentium Intel Core2 Quad 2.66 GHz processor. Since our 
aim is to estimate the scaling of quantities right at the critical point, simulations must be 
performed at the critical temperature of the model. We used the estimates for the critical 
temperature given in [TO]. In terms of (3 = 1/kT, where k is the Boltzmann constant, 
these are pc = 0.149695, /3c = 0.188864 and /3c = 0.317368, for p = 1, p = 0.8 and p = 0.5, 
respectively. 

In addition we have simulated the dilution p = 0.650 at Pc = 0.235049 [lU] with 
the same statistics used for the other dilutions but we have found that the behaviour 
of the observables differs from that expected. In Fig. [2] is shown the deviation of the 
leading scaling behaviour of the susceptibility. We have rechecked this point starting 
from different initial configurations and even using different random number generators. 
This is probably due to a biased estimation of the critical temperature in [T^ . For this 
reason we omit p = 0.650 from our analysis. 
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Tabic 3: Simulation details for each spin concentration p and system size L. Here, A^voifr 
denotes the number of Wolff updates between consecutive measures, A^d is the number of 
dropped measurements at the beginning of the Monte Carlo history (in units of 10^), and 
A^m is the total number of measurements performed after thermalization in units of 10^. 
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Figure 2: (Color online) Comparative behaviour of the susceptibility for p = 0.650 at 
Pc = 0.235049 [in]. Note the strong deviation at this dilution from the expected leading 
behaviour [x ~ -^^)- The point size is in every case bigger than the corresponding error 
bar. 



5 Analysis of the Numerical Results 

To establish confidence in the present approach, the pure system is analysed first to 
test whether the method employed successfully quantitatively identifies the logarithmic 
corrections which are well established there. 



5.1 The Pure Case p = 1 

The scaling and FSS of the pure model {p = 
heat FSS behaviour is given by 



1) are well understood [T0l[T2]. The specific 



Cl(0) ~ (hiL)" ~ (lnL)i/3^ (5.1) 

up to additive corrections. Fitting to this form for a over the full data set 8 < L < 48, 
one finds the estimate a = 0.42(4) with a goodness of fit given by a x^/d.o.f. (chi-squared 
per degree of freedom) of 1.3. The estimate is two standard deviations away of the known 
value 1/3. As elsewhere in this analysis, inclusion of sub- leading scaling correction terms 
in the fits does not ameliorate this result, which is similar to that reported in [TO] . 

The FSS for the susceptibility is given in (13. 7p with ( = 1/2. Fitting to the leading 
form 

Xl(0)~L^ (5.2) 
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gives 7/z/ = 2.16(1) for 8 < L < 48 and 7/z/ = 2.13(2) for 12 < L < 32, the difference 
from the theoretical value 7/z/ = 2 being ascribable to the presence of the logarithmic 
correction term. Accepting this mean- field value for 7/z/ and fitting to 

XL{0)^L\\nL)^, (5.3) 

gives the estimate ( = 0.48 ± 0.02 in the range S < L < 48, albeit with x^/d.o.f. = 4.1 

The FSS for the individual Lee- Yang zeros is given in (13.91) with p = — 1/4 in the pure 
case. Fitting to the leading form 

rj{L) ~ L-^ (5.4) 

gives A/z/ = 3.074(5) for 8 < L < 48, the difference from the theoretical mean-field value 
A/z/ = 3 being due to the corrections. Accepting this value and fitting to 

rj(L) ~ L-=^(lnL)^ (5.5) 

gives p = —0.22(2) in the range 8 < L < 48. This estimate is compatible with the known 
value p = —1/4. As one would expect, the higher zeros yield less accurate estimates (as 
they are further from the real simulation points) with p = —0.18(3) , p = —0.17(7) and 
p = —0.10(14) from the second, third and fourth zeros respectively. These estimates are 
listed in Table HI 

Having established that the numerics gives reasonable agreement with the pure theory 
at the leading and the logarithmic levels, we now perform a similar analysis in the presence 
of disorder. 

5.2 The Diluted Cases p = 0.8 and p = 0.5 

Since the Ansatz (13.61) for the specific heat in the disordered systems is rather more com- 
plex than that of the pure case (15. ip . we begin the p 7^ 1 analyses with the susceptibility 
and the Lee- Yang zeros. It will turn out that our analyses will reinforce the analytical 
predictions that scaling is governed by the Gaussian fixed point and that the logarith- 
mic corrections in the RSIM differ from those in the pure model. Indeed, the results for 
the zeros will be seen to be broadly compatible with the analytic predictions contained 

in iiiTiiHiEino]. 

For the weaker dilution value p = 0.8, a fit using all lattice sizes to the leading form 
(15.21) for the susceptibility yields the estimate 7/z/ = 2.14±0.01. 

Ascribing the difference from the Gaussian value 7/z/ = 2 as being due to the correction 
terms and, as in the pure case, and fitting to (15. 3p . one finds an estimate for the correction 
exponent ( = 0.39(3) for 8 < L < 48. This values is between the pure value ( = 0.5 and 
the theoretical estimates for the diluted value which give ( ~ 0.25 to 0.26. Thus, while 
the FSS for the susceptibility does not capture the theoretical estimates for the diluted 
case, the fitted values have moved away from the pure value and towards the lower value 
listed in Table [2l As elsewhere in this work, the inclusion of scaling corrections does not 
alter these results significantly. 

A similar analysis for the FSS of the susceptibility at the stronger dilution value p = 0.5 
gives similar results: the leading form (15. 2p yields an estimate 7/z/ = 2.13 ± 0.02 with 
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Table 4: FSS estimates for the various dilution values, using a range of lattice sizes. 
The susceptibility is expected to scale as xl ~ L^(lnL)'' and the Lee- Yang zeros as 
Tj ~ L^^(lnL)^, where ( ~ 0.25 to 0.259 and p ^ —0.125 to —0.130. (For comparison, 
the pure theory with p = 1 has ( = 1/2 and p = —1/4.) 
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0.42(4) 
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-0.17(5) 


-0.18(4) 
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L = 8 - 48 


0.37(4) 


-0.20(4) 


-0.22(4) 


-0.21(4) 


-0.21(4) 
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-0.20(5) 


-0.18(5) 


-0.19(5) 
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In L 



1 1.2 
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Figure 3: (Color online) (a): FSS plot for at p = 0.8 (circles) and p = 0.5 (triangles) 
at the critical point. The slopes of the fitted solid and dashed lines give estimates for 7/1/ 
of 2.14(1) and 2.13(2), respectively, (b): Plot of In^L - 21nL against In (InL) at p = 0.8 
(circles) and p = 0.5 (triangles) giving slopes 0.39(3) and 0.37(4), respectively, indicating 
slow crossover of multiplicative logarithmic corrections from the pure case (where ( = 0.5) 
to the diluted case, where the theoretical value is C ~ 0.13. 



a goodness of fit given by x^/d.o.f. = 0.4. Ascribing the difference from the mean-field 
value 7/z/ = 2 as being due to the logarithmic corrections, and fitting to (15.31) . one obtains 
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the estimate ( = 0.37(4) for 8 < L < 48. Again this result is between the theoretical 
predictions for the pure {( = 0.5) and diluted {( ~ 0.25 to 0.26) cases. These results 
are summarised in Table HI together with results obtained from the same fits with the 
smallest lattices removed. The best fits of the susceptibility can be seen in Fig. [31 

Since in each of the diluted cases, the results for susceptibility lie between what is 
expected for the pure theory and for the diluted theories, we appeal to the Lee- Yang 
zeros, as our collective experience indicates that they give a cleaner signal. 




U I I I I I I I I I I I I I I I I I I I U \J ^ t I I I I I I I I I I I I l_=l 

2 2.5 3 3.5 4 ' 0.8 1 1.2 1.4 
In L In (In L) 

Figure 4: (Color online) (a): FSS plot for the Yang-Lee edge at p = 0.8 (circles) and 
p = 0.5 (triangles). The slopes of the fitted solid and dashed lines give estimates for A/z/ 
of 3.055(4) and 3.07(2), respectively, (b): Plot of Inri + 3 InL against In (InL) at p = 0.8 
(circles) and p = 0.5 (triangles). Fits in the range L = 12 to L = 48 (plotted) give slopes 
—0.17(4) and —0.16(5), compatible with the predictions ranging from p ^ —0.125 and 
p ^ —0.13 in the literature. (For comparison, in the pure model, p = —1/4). 



The leading behaviour is firstly examined by fitting each of the first four Lee- Yang 
zeros to (15.41) . For the weaker dilution given hj p = 0.8, one obtains A/z/ = 3.055(8), 
3,056(9), 3.069(11) and 3.060(10) from fits to the first, second, third and fourth zeros 
respectively, using all lattice sizes. The equivalent results for the stronger dilution value 
p = 0.5 are A/z/ = 3.068(13), 3, 071(15), 3.072(12) and 3.071(11), respectively. All fits are 
of good quality with acceptable values of x^/d.o.f ., which we refrain from detailing. Again, 
these are interpreted as being supportive of the mean- field leading behaviour 7/z/ = 3 with 
logarithmic corrections. 

The logarithmic-correction exponents are estimated by fitting to 05.51) . with the various 
theories indicating that p = —0.125 to —0.13. The strongest evidence supporting this 
comes, as it should, from the first zero (the Yang-Lee edge) for p = 0.8, which yields 
the estimate p = —0.15(2) (with x^/d.o.f. = 0.6). As in the pure case, and as expected, 
estimates for p deteriorate as higher-index zeros are used. Dropping the smallest lattices 
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from the analysis, however, leads to these estimates for p more compatible with [6l [71 [8], 
m [To] . These results are summarised in Table [H 

The equivalent analysis for the stronger dilution value p = 0.5 is less clear, with an 
estimate p = —0.20(4) coming from the first zero when all lattices are included in the fit 
(with x^/d.o.f. = 1.1). Dropping the smallest lattices, however, gives p = —0.16(5) (with 
X^/d.o.f. = 0.9), closer to the values coming from [6], [71 [H [9l [10]. Similar results are 
obtained for the higher zeros and these are also summarised in Table [H The best fits for 
the first zero are shown in Fig. [H 

As a final check of the reliability of our results we have used the spectral energy 
method [Tl| to reweight the data obtained at Pc to Pc ± ^Pc (taken again from (TU]) 
obtaining that the new data sets are fully supportive of the previous result^. 

Having checked that the leading FSS behaviour corresponds to that originating in the 
Gaussian fixed point, and that the logarithmic corrections to scaling are different to those 
in the pure model, and moreover (at least in the case of the Yang-Lee edge) are broadly 
compatible with the predictions from the literature [SI [TJ [HI [H [ID], we now attempt to 
distinguish between these broad predictions. To this end we turn to the specific heat. 




10- 



10 20 ^ 30 40 50 

Figure 5: (Color online) The specific heat for p = 0.8 (circles) and p = 0.5 (triangles). 
The error bars are in every case smaller than the point size. The solid and dashed curves 
are best fits to the Ansatz (13.61) . with a = 1/2. 

Having established confidence that the mean-field values 7=1 and A = 3/2 hold in 

^ When doing /3-extrapolations in a disordered model one should be careful and take into account the 
bias induced by the finite number of measures (see the discussions in [lOl [15] ) . We have followed the 
recipe provided in '10] to perform the extrapolation to infinite number of measures per sample. 
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the 4D RSIM, we may use the scahng relation a = 2 — 2A + 7 to estabhsh the mean- 
field value a = too. The Ansatz (13. 6p for the specific heat may now be used. This 
contains information which can be used to discriminate between some of the scenarios in 
the literature. From Table [H there is a striking difference between the estimates for the 
specific- heat logarithmic-correction exponent a coming from [71 [H] and from [HI El 110] • 
While the former have relatively large values of a, the latter agree on d = 0.5. The 
simulated values of the specific heat at p = 0.8 and p = 0.5 are plotted in Fig. [51 The 
slope of the full specific-heat curve (13. 6p is 




^ = [A - C.m^r^ I 1 - ^Sf=^ 1 • (5-6) 

This vanishes when Cl{0) = A and when y/hiL = a ^53/ 12. The first of these is the 
asymptote L — 00, from which A can be determined for each dilution. The second 
occurrence of zero slope is for quite small lattice sizes, i.e., beneath lattice size L = 8. 



Therefore a ^ y 53/12vin8 ~ 0.7, excluding the values a ~ 1.237 and a ~ 1.246 given 
in [3 [9]. In fact, a best fit to the Ansatz ([3lD gives A = 49(11), B' = 66(22) and 
a = 0.46(18) for p = 0.8 and A = 10(5), B' = 7(10) and a = 0.7(3) for p = 0.5 with 
L > 8. Fixing a = 1/2 in each case gives A = 52(2), B' = 72(5) for p = 0.8 and 
A = 9(3), B' = 5(1) for p = 0.5, and these curves are plotted along with the specific heat 
measurements in Fig. [5l Fixing the correction exponent a to the value given in [3, [9], on 
the other hand, yields a best-fit value of B' which is negative in each case, contradictory 
to [H H] . Thus we can deem these values to be unlikely. 



6 Conclusions 

Numerical measurements of the leading critical exponents in the 4D RSIM are presented, 
confirming that the phase transition in this model is governed by the Gaussian fixed 
point. We then turn to the corrections to scaling, for which there exist five distinct sets of 
predictions in the literature [SI 13, [SI [SI [ID] • The scaling relations for logarithmic corrections 
are used to render complete these sets and their counterparts for finite-size systems are 
given. 

The measured values of the susceptibility FSS correction exponent ( for the site-diluted 
model, lie between the known value for the pure model and the theoretical estimates 
coming from P [3 [SI [SI [ID] for the disordered system. While this result illustrates slow 
crossover of the susceptibility, the lowest lying Lee- Yang zeros give a cleaner signal and 
the measured value for their logarithmic correction exponents are indeed compatible with 
the theories. 

To discriminate between the five theories, the detailed finite-size scaling behaviour 
of the specific heat is also examined. The analysis is clearly in favour of the analytical 
predictions of [HI [SI HD] over those of [71 [S]- This is contrary to expectation as the for- 
mer involve only two loops in the perturbative RG expansion, while the latter take the 
expansion to three loops in the beta function. 
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